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1 PRESENTACIÓN 


El presente proyecto se sustenta en el empleo de herramientas de Teledetección para 
identificar los componentes ambientales más relevantes en la caracterización del paisaje natural 
en el cual se desarrolla la ruta QHAPAQ ÑAN. 

En términos generales, el concepto de Teledetección puede ser entendido como el estudio o el 
análisis de cualquier objeto sin entrar en contacto con él. Al respecto la NASA 1 define como 
Teledetección: 


“La adquisición y registro de datos o información de alguna o varias características 
de un fenómeno, objeto o material, a través de dispositivos que no entran en 
contacto con ellos. ” 


Otro de los elementos que caracterizan la Teledetección es el empleo de diversas plataformas, 
ya sean terrestres, marítimas, áreas o espaciales, las cuales entregan diversas gamas de 
aplicaciones. Al respecto, el proyecto se estructura a partir del uso de imágenes obtenidas a 
través de plataformas satelitales, específicamente el Sensor Koreano Kompsat II. 

La utilización de la Teledetección en el estudio y análisis del territorio presenta un sin número 
de ventajas, entre ellas destacan: la relativa rapidez y bajo costo de mantener actualizada 
información de amplias áreas geográficas, permite obtener información de sectores 
inaccesibles y principalmente la información puede ser integrada con otras bases de datos 
georreferenciadas a través de un Sistema de Información Geográfica. 

En este contexto, el Proyecto QHAPAQ ÑAN cuenta con información satelital, cuyas 
características técnicas (resolución espectral y espacial), permiten en líneas generales una 
caracterización adecuada de los elementos abióticos, bióticos y antrópicos existentes en la ruta 
patrimonial. 

El siguiente informe, presenta los alcances del proyecto “Levantamientos terrestres para los 
tramos del itinerario y paisaje cultural del Qhapaq Ñan-Chile, análisis y productos obtenidos de 
la teledetección”, cuyo desarrollo fue llevado a cabo por Andalué Ambiental S.A. 

El encargo realizado busca generar información digital a escala 1:5000 de diversos elementos 
en 13 sitios seleccionados por la contraparte, los cuales forman parte de los distintos itinerarios 
de Qhapaq Ñan y área de Paisaje Cultural a saber: 


1 Fuente NASA sitio web http://rst.gsfc.nasa.gov/lntro/Part2_1 .html 
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Tramo o subtramo 

Sitio 

Subtramo 1: Putre-Zapahuira 

Sitio N S 1: Que abarca desde Putre hasta 
Zapahuira 

Tramo 2: Alto Loa. Miño - Lasaña 

Sitio N 2 2: al sur de Embalse Conchi 

Sitio N 2 3: al sur de Lasaña 

Áreas de Paisaje Cultural: Ayquina, Toconce, 
Caspana, Turi, Cupo - Paniri 

Sitio N 2 4: Abarca la totalidad de esta área 

Subtramo 3:Cupo-lncahuasi; Catarpe 

No tiene sitios asociados 

Subtramo 4A: Despoblado de Atacama: Camar - 
Peine 

Sitio N 2 5: sector de Camar 

Sitio N 2 6: sector de Socaire, Mal paso 

Sitio N 2 7: sector de Peine 

Subramos 4B: Despoblado de Atacama: Portal 
del Inca-Sierra de Capis 

Sitio N 2 8: sector Portal del INca 

Sitio N 2 9: sector Finca Chañaral 

Sitio N 2 10: sector pueblo Inca de Oro 

Subtramo 5: Punta Brava-Viña del Cerro; Iglesia 
Colorada 

Sitio N 2 11: sector Punta Brava 

Sitio N 2 12: sector Las Puertas 

Sitio N 2 13: area de Viña del cerro 


En estos sitios se trabajó con imágenes KOMPSAT-2, sobre las cuales se desarrollaron 
diversos procesos tendientes a corregir distorsiones propias de la toma de datos 


El trabajo, fue desarrollado en dos etapas: 

• La primera considera las correcciones necesarias para comenzar el procesamiento de 
las imágenes satelitales, cuyos productos correspondieron a los mosaicos de imágenes, 
las imágenes seucolor, NVDI, reflectancia y falso color convencional. 

• La segunda etapa se constituyó de los productos finales del procesamiento de imágenes 
satelitales y la correspondiente generación de composiciones cartográficas. 

Los productos anteriores, además de la resolución de las observaciones realizadas por parte de 
la contraparte técnica, forman parte integrante de este documento final. 
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2 CARACTERÍSTICAS DE LAS IMÁGENES KOMPSAT-2 II. 

El siguiente proyecto se ha desarrollado sobre la base de imágenes KOMPSAT-2, satélite 
Coreano puesto en órbita el 28 de julio de 2006 2 . Las características de resolución espectral, 
espacial y sus aplicaciones se sintetizan en la siguiente tabla: 


Tabla 1. Características de las imágenes KOMPSAT-2 


Sensor 

Rango 

espectral 

de 

la banda 

Banda 

Ancho de 
la 

banda (pm) 

Resolución 

espacial 

(m) 

Utilidades 

Nombre 

Tipo 

KOMPSAT-2 

Pancromático 

Visible - 

infrarrojo 

cercano 

Banda P 

0,45-0,9 

1 

Localizan, identifica y 
miden accidentes 

superficiales y objetos, 
principalmente por su 
apariencia física 

Multiespectral 

Azul 

Banda 1 

0,45 - 0,52 

4 

Cartografía de aguas 
someras, diferenciación 
de suelo y vegetación 

Verde 

Banda 2 

0,51 -0,6 

4 

Diferenciación de 

vegetación por su salud 

Rojo 

Banda 3 

0,63-0,7 

4 

Diferenciación de 

vegetación por su 

especie? 

Infrarrojo 

Cercano 

Banda 4 

0,76-0,85 

4 

Cartografía de la 

vegetación 

Cartografía del vigor/ 
salud de la vegetación, 
diferenciación de 

vegetación por especies? 


Fuente: Tomado desde bases de licitación elaboradas por Ministerio de Interior, Subsecretaría de Desarrollo Regional y 
Administrativo 


En la tabla 2 se muestran el ancho del espectro y la constante o irradiancia solar para cada una 
de las bandas multiespectrales. 


2 Lee, D., 2008, Image Restoration of Calibration and Validation for KOMPSAT-2, ISPRS2008 
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Tabla 2. Características Espectrales del Sensor KOMPSAT II 


Banda 

Rango Espectral (pm) 

Constante Solar (W/m 2 /ster/pm) 

1 

0,45 - 0,52 

1918,6 

2 

0,51 -0,6 

1840,4 

3 

0,63-0,7 

1545,0 

4 

0,76-0,85 

1073,5 


Nota: pm: micrometro; W: Watts; ster: Stereo Radian, 
Fuente: R. Santer, 2007 3 


3 R. Santer, 2007. CalVal Support System, Brockman Consultant. 
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Las imágenes abarcan los años 2008 al 2009 En la tabla 3 se presenta un resumen de las 
principales características de las 26 escenas proporcionadas. 


Tabla 3. Escenas Kompsat II Procesadas. 


Sitio 

Fecha 

Coordenada Central 

Hora UTC 

Distancia 
Tierra - Sol 

Latitud 

Longitud 

Sltio 01a 

09-11-2008 

-18,3810 

-69,5903 

15,2954 

0,9902 

Sltio 01b 

09-11-2008 

-18,3810 

-69,5903 

15,2954 

0,9902 

Sltlo 01c 

09-11-2008 

-18,1997 

-69,5778 

15,2954 

0,9902 

Sltlo 02a 

25-07-2008 

-22,0591 

-68,6280 

15,3741 

1,0157 

Sltlo 02b 

25-07-2008 

-22,1563 

-68,6341 

15,3741 

1,0157 

Sltlo 03a 

09-11-2008 

-22,3035 

-68,6950 

15,2954 

0,9902 

Sltlo 03b 

25-07-2008 

-22,2633 

-68,6091 

15,3741 

1,0157 

Sltlo 03c 

25-07-2008 

-22,2633 

-68,6091 

15,3741 

1,0157 

Sltlo 04a 

17-10-2008 

-22,1295 

-68,3615 

15,2647 

0,9962 

Sltlo 04b 

17-10-2008 

-22,2633 

-68,3307 

15,2647 

0,9962 

Sltlo 04c 

17-10-2008 

-22,3971 

-68,2998 

15,2647 

0,9962 

Sltio 04e 

19-05-2009 

-22,1390 

-68,2826 

15,2647 

1,0116 

Sltlo 04f 

19-05-2009 

-22,2383 

-68,2372 

15,2647 

1,0116 

Sltlo 04g 

27-05-2009 

-22,3361 

-68,2233 

15,3741 

1,0132 

Sltio 04h 

17-01-2008 

-22,2724 

-68,1121 

15,3741 

0,9838 

Sltlo 05 

06-06-2009 

-23,4099 

-67,9799 

15,3741 

1,0147 

Slto 06 

06-06-2009 

-23,5791 

-67,9076 

15,3741 

1,0147 

Sltio 07a 

17-10-2008 

-23,6009 

-68,0208 

15,2339 

0,9962 

Sltlo 07b 

17-10-2008 

-23,7347 

-67,9896 

15,2339 

0,9962 

Sltio 07c 

17-10-2008 

-23,6794 

-68,0755 

15,2339 

0,9962 

Sitio 08 

17-11-2008 

-26,2742 

-69,5768 

15,4675 

0,9885 

Sltlo 09 

12-05-2009 

-26,6482 

-69,9052 

15,4306 

1,0101 

Sltio 10 

12-05-2009 

-26,7818 

-69,8699 

15,4400 

1,0101 

Sltlo 11 

26-02-2008 

-22,3035 

-68,6950 

15,6025 

0,9900 

Sltio 12 

11-09-2007 

-27,8635 

-70,0520 

15,7000 

1,0064 

Sltio_13 

11-09-2007 

-27,8925 

-70,0427 

15,7010 

1,0064 


Fuente: Elaboración propia. 
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2.1 Imágenes Base Para la Entrega de Resultados 

La generación de imágenes de Reflectancia, NDVI, Seudo Color (Bandas 321) y Falso Color 
(Bandas 432), se desarrollaron a partir de las imágenes proporcionadas por el CMN, en la tabla 
4 se detalla la imagen base para cada resultado. 


Tabla 4. Archivos e Imágenes Base y Resultado. 


Archivo e Imagen Base 

Carpeta e Imagen Resultado 

Archivo o 
Carpeta 

Imagen Base 

Carpeta de 
Resultado 

Imagen Resultado 

Archivo 
Imágenes.rar 

Carpeta DEM 

Imágenes Kompsat II 
Formato Geotlff Sitio_01 a 
Sitio_13 bandas 4, 3, 2 y 1 
con resolución espacial de 

4 metro. 

Imágenes DEM (altitudes) 
en formato Geotiff con una 
resolución espacial de 10 
metros. (Utilizadas en la 
corrección topográfica 
para la generación de las 
Imágenes 

Reflectancia 

Imágenes Formato Idrisl para 
cada una de las cuatro bandas 
en los 13 sitios. 

NDVI 

Imágenes Formato Idrisl para 
cada una de las cuatro bandas 
en los 13 sitios. 

Falso Color 
Convencional 

Imágenes Formato Idrisl para 
cada una de las cuatro bandas 
en los 13 sitios. 

Archivo 

Strech.rar 

Imágenes Formato Geotiff 
Sitlo_01 a Sitio_13 bandas 
3, 2 y 1 con resolución 
espacial de 1 metro. 

Pansharp 

Imágenes Formato Idrisi 

Sitlo_01 a Sltio_13 bandas 3, 2 
y 1 con resolución espacial de 

1 metro. 
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3 METODOLOGÍA 


3.1 Correcciones Atmosféricas, Radiométricas y Topográficas. 

Estos procedimientos buscan estandarizar la información contenida en todas las imágenes, 
siendo un paso necesario para poder efectuar una comparación entre bandas para las 
imágenes procesadas. Al final de este procedimiento cada una de las cuatro bandas 
procesadas, presentará valores en un rango de 0 a 100, lo que expresará la cantidad de 
energía reflejada por la superficie terrestre al momento de la captura de la imagen satelital. 


3.1.1 Correcciones A tmosféricas 

Para la corrección atmosférica de cada una de las cuatro bandas procesadas, se aplicó el 
método de Sustracción del Objeto Oscuro (Song, et al. 2001), éste procedimiento ampliamente 
utilizado para clasificaciones de imágenes, consiste en determinar a través del histograma de la 
escena su valor mínimo ya sea de nivel digital (nd) o radiancia, ese valor luego es sustraído a la 
escena completa (Chávez, 1989), este proceso se realiza en cada una de las bandas, siendo 
mayores las correcciones en las bandas azul y verde, disminuyendo luego en las bandas rojo e 
infrarrojo (ver ejemplo Figuras 1 y 2). 


Figura 1. Ejemplo de Histograma Imagen Sitio 02a Banda 1 (Azul) 

Histogram of 'sitio_02a_band1 ’ using 'sitio_02a_masc' as mask 

Summary Statistics 

Classwidth Mean Actual min Actual max N 

1 j 318.808 _|_0_ 510 1070489 

14000 
12000 
10000 
8000 
6000 
4000 
2000 
0 

29 59 89 119 149 179 209 239 269 299 329 359 389 419 449 479 509 


Fuente: Elaboración propia. 



Std deviation 
40.444 
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Figura 2. Histograma Imagen Sitio 02a Banda 1 (Azul) luego de Aplicar el Método del Pixel Más 
Oscuro. 


Histogiam of 'sltio_02a_band1 af using 'sitio_02a_masc' as mask 

Summary Statistics 




Fuente: Elaboración propia. 


Tal como se aprecia en las Figuras 1 y 2, se presenta el histograma de la imagen sitio_02a en 
su banda 1, en el eje de las abscisas se presentan los niveles digitales y en el eje de las 
ordenadas la cantidad de celdas. De este modo la Figura 1, muestra la imagen con su 
distribución de celdas para cada nivel digital en su estado original, para luego en la Figura 2 
indicar como se reducen los valores de los niveles digitales producto de la sustracción del pixel 
más oscuro, el cual representa la influencia de las condiciones atmosféricas en la captura de la 
imagen satelital. 


3.1.2 Correcciones Radiométricas 


La calibración radiométrica consiste en la transformación entre los valores correspondientes al 
nivel digital de cada pixel (información de como el sensor almacena los datos recibidos) a 
parámetros físicos (valores referidos a la reflectancia) de la cubierta terrestre. 


El cálculo de reflectancia se efectuó implementando el algoritmo COST de Chávez (1996), el 
cual se expresa: 

o wat- x(rBandN)xD 2 

P ütl EBandN x (COS ((90 - 0) x n /180)) x r 


Andalué 


INFORME FINAL: “Levantamientos terrestres para los tramos del itinerario y paisaje cultural del 
Qhapaq Ñan-Chile, análisis y productos obtenidos de la teledetección" 













CONSEJO DE MONUMENTOS NACIONALES 


ANDALUE AMBIENTAL S.A. 


Donde: 

pBandN . R e f| ec tar,c¡a para la banda N. 

rBandN ■ Radiancia con corrección atmosférica para la banda N. 

D : Distancia Tierra-Sol para el día de la imagen. 

EBandN : irradiancia solar para la banda N. 

^ : Ángulo de elevación solar para la hora y día de la imagen. 

T : Transmitancia atmosférica expresada como (C0S'((9O-6 , )x;r/18O)) 

Cabe señalar que los valores para EBandN D y # se muestran en las tablas 1 y 2 
respectivamente. 


Figura 3a. Imagen Sitio 03a Banda 1 Figura 3b. Imagen Sitio 03a Banda 1 

Previo al Cálculo de Reflectancia Luego del Cálculo de Reflectancia 



Nota: Los valores de reflectancia se expresan en el rango de 0 a 100. 
Fuente: Elaboración propia. 
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Las Figuras 3a y 3b muestran un ejemplo de la imagen del sitio 03a para su banda 1 (espectro 
azul) luego de realizada la corrección atmosférica y posteriormente de aplicar el algoritmo de 
calculo de reflectancia. De este modo se aprecia como en la Figura 3a los valores de niveles 
digitales fluctuaban entre -150 y 844, para posteriormente ser transformados a valores físicos 
de reflectancia en el rango de 0 a 41,7, lo cual permite comparar este parámetro para las otras 
tres bandas (verde, rojo e infrarrojo cercano) y poder construir la firma o curva espectral de los 
objetos que serán clasificados en la etapa 2 del presente estudio. 


3.1.3 Correcciones Topográficas 

Para minimizar la influencia del relieve en los valores de reflectancia registrados en cada una de 
la escena, se realizó una corrección topográfica, donde a partir del Modelo Elevación Digital 
(DEM), el día, la hora y la ubicación geográfica se genera un hillshade o cobertura de sombra 
para cada escena. 

El hillshade se generó con el módulo SURFACE del programa Idrisi Andes, ingresando el DEM 
de 10 metros de pixel, junto con el ángulo de elevación solar y el ángulo azimutal solar, valores 
que son calculados a partir de la hora UTC y la coordenada central de la escena respectiva. El 
resultado es una imagen de sombras con valores de 0 a 1, indicando los sectores de menor a 
mayor iluminación, tal como se muestra en las figuras 4a y 4b. 
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Figura 4a 

DEM Sector Inkahuasi 


Figura 4b 
Hillshade 



Fuente: Elaboración propia. 
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De esta manera, se aplico el procedimiento de Calibración Topográfica de Idrisi Andes, cuyo 
resultado es una estimación de la reflectancia que captaría el sensor si el terreno fuese 
completamente plano, el procedimiento contempló los siguientes pasos: 


1. A través del programa HISTO se calculó el valor medio de reflectancia para cada banda 
a ser corregida. 

2. Usando el DEM para la imagen del área, se creó un hillshade analítico por medio del 
Módulo SURFACE, tal como el presentando en la figura 4b. 

3. Con el módulo CALIBRATE, se procedió a realizar una regresión entre la imagen 
hillshade y la imagen de la banda a ser corregida. En este punto lo que se persigue es 
establecer que tendencia se presenta entre los valores de reflectancia de la banda 
analizada y las condiciones de sombra que presenta el relieve para el sitio analizado. 
Esto generó una nueva imagen la que se utilizó en el siguiente paso. 

4. A través del módulo Image Calculator, se procedió a realizar una resta entre la imagen 
de la banda a corregir y la imagen generada en el punto anterior, y luego se sumo el 
valor promedio calculado en el punto 1. De este modo se obtiene una aproximación a las 
condiciones reales de reflectancia que registró el sensor satelital al momento de captura 
la respectiva escena. 


Un ejemplo del resultado del procedimiento anterior se muestra en las Figuras 5a y 5b, donde 
se puede apreciar que los valores de reflectancia presentaron un ligero cambio, lo cual se 
explica porque las condiciones topográficas del sitio en cuestión son bastante planas tal como 
se constato en la Figura 4b, donde sólo en una quebrada al centro de la imagen se presentan 
condiciones de sombras que pueden llegar a reducir la iluminación hasta en un 30%, ya que la 
imagen presenta valores entre 0,68 y 1. En consecuencia en aquellos sectores de relieve más 
rugoso o montañoso, las diferencias de reflectancia entre una imagen corregida y no corregida 
topográficamente son indudablemente más significativas. 
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Figura 5a 

Imagen Sitio 03a Banda 1 
Previo a la Corrección Topográfica 


Figura 5b 

Imagen Sitio 03a Banda 1 
Luego de la Corrección Topográfica 



Fuente: Elaboración propia. 


A modo de síntesis e ilustrando de una manera simplificada los procesos anteriormente 
expuestos, se presenta en la Figura 6 un esquema cartográfico con los diversos procesos e 
imágenes intermedias generadas, para obtener el producto de imagen de reflectancia para cada 
una de las cuatro bandas (azul, verde, roja e infrarrojo cercano) en los sitios analizados. 
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Figura 6: Esquema Cartográfico (Algoritmo de obtención de imagen corregida) 
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3.2 Generación de índices de Vegetación 


El proceso de generación de índices de Vegetación tiene por objetivo detectar zonas de 
vegetación y establecer cuantitativamente la vigorosidad de ésta, así como develar las zonas de 
ausencia de cobertura vegetal. En primer lugar se aplicó el algoritmo del índice Normalizado de 
Vegetación (NDVI en su sigla en inglés), ampliamente utilizado en el estudio de la vegetación 
(Lillesand y Kiefer, 2004). El NDVI busca registrar la diferencia espectral existente entre la 
banda del infrarrojo cercano y el rojo, donde el resultado fluctúa entre -1 y 1, siendo los valores 
mayores a cero los que indican la presencia creciente de vegetación. El algoritmo se 
implemento de la siguiente forma: 


BandalRC - BandaR 


NDVI = 


BandalRC + BandaR 


La implementación del NDVI, permitirá desarrollar dos lineas de análisis. La primera centrada 
en las areas de vegetación, donde posteriormente se discriminará entre la existencia de 
cultivos, vegas, bofedales y otras cubiertas vegetales. En tanto, en paralelo se podrá analizar 
todas aquellas áreas que no corresponden a coberturas vegetacionales, vale decir: nubes, 
rocas, suelos desnudos, nieves, ríos, lagos, caminos, edificaciones entre otros. Como la 
vegetación del área de estudio suele ser muy escasa en comparación al suelo desnudo (por 
tratarse por una zona árida) utilizamos el Indice de Vegetación Ajustado al Suelo (SAVI). Este 
se obtuvo a partir de la siguiente ecuación (Huete, 1988) 


(.BandalRC - BandaR)(\ + L) 
{BandalRC + BandaR + L ) 


Siendo L el factor de correción utilizado para disminuir el efecto de la reflectancia del suelo 
desnudo, usándose comunmente un valor de 0,5 (Huete, 1988). En esta oportunidad, ocupamos 
un factor L de 1.0 por tratarse de zonas áridas. 


Ambos procedimientos se encuentran en el módulo VEGINDEX de IDRISI Andes. 
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Figura 6a Figura 6b 

Imagen Sitio 07a Imagen Sitio 07a 

NDVI SAVI 



En síntesis la Figura 6a y 6b muestran los resultados al implementar los algoritmos NDVI y 
SAVI para el caso del Sitio 07. En la Figura 6a se aprecia una mayor diversidad de rango y 
donde la vegetación detectada al sur de la imagen, alcanzan los valores máximos de NDVI (1), 
no obstante estos valores pueden estar influenciados por la aridez del terreno, sobre estimando 
su valor. En contrapartida la Figura 6b muestra que al aplicar el algoritmo SAVI, se genera una 
imagen con un rango más reducido de valores y donde el valor máximo alcanza a 0.29, lo que 
evidentemente atenúa la influencia del suelo desnudo, condición que prevalece en zonas áridas 
como es el caso del área de estudio del presente proyecto. 
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3.3 Creación de Falso y Seudo Colores 

Con el módulo COMPOSITE de IDRISI Andes, hemos generado imágenes falso color 
convencional o 432. En ella la Banda 4 (Infrarrojo cercano) ocupa el cañón Rojo, la banda 3 
(Rojo) el cañón Verde y la Banda 2 (Verde) en el cañón Azul del sistema de colores RGB. 
Todas estas escenas poseen una resolución espacial de 4 metros, y se utilizarán principalmente 
para la detección de zonas de vegetación, que aparecen en rojo. 

Con el módulo PANSHARP de Idrisi Andes, hemos generado un mosaico de imágenes a mayor 
resolución (1 m) a partir de utilizar una de las bandas RGB de las imágenes del archivo Strech 
proporcionado por CMN las cuales se presentaban con pixel de 1 metro 4 . Las bandas azul, 
verde y roja ocupan cada una su cañón respectivo del sistema RGB, dando como resultado una 
imagen de seudo color natural. Con ella realizaremos inspección en detalle de los elementos 
geográficos presentes en la superficie del área de estudio. 



Figura 7a 
Imagen Sitio 11 
Falso Color 432 


Figura 7b 
Imagen Sitio 11 
Seudo Color 321 


La creación de estos Falsos o Seudo Colores lo que persiguen es facilitar la interpretación 
visual de las imágenes en las etapas posteriores de este estudio. Por un lado el falso color 
convencional resalta en tonos rojos la vegetación, siendo este un color fácilmente reconocible 


4 Cabe destacar que no se tuvo acceso a las imágenes pancromáticas ortorrectificadas para cada uno de los sitios, 
por ende el procedimiento de Pansharp fue realizado solo con las imágenes strech. 


% 
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por el ser humano y por ende ayuda en la clasificación de objetos y coberturas vegetacionales 
presentes en la imagen. En tanto, el seudo color o también llamado falso color real o natural, 
simula los tonos que observa normalmente una persona a partir de las condiciones que 
presenta el territorio, tal como sería el caso si observase en pleno vuelo en un avión. 

Ambas técnicas permiten identificar y discriminar visualmente las características de los 
elementos presentes en la imagen tales como: rocas, suelos desnudos, cuerpos de agua, 
cultivos entre otros, tarea especialmente sensible al momento de tomar muestra o sitios de 
entrenamiento para clasificar la imagen y poder generar los mapas temáticos respectivos. 


4 SIGNATURAS ESPECTRALES E INTERPRETACION DE IMAGENES 

El primer paso para la clasificación temática de las 26 escenas fue el análisis visual de las 
imágenes de seudo color (imágenes compuestas con las bandas 3-2-1) y falso color 
convencional (imágenes compuestas con las bandas 4-3-2) para realizar la digitalización de 
sitios de muestreo o “entrenamiento” (figura 8). 


Figura 8 Sitios de entrenamiento para la escena del sitio 12 (Los Loros) 
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A partir de los sitios de muéstreos se obtuvo la firma espectral de los distintos tipos de 
coberturas (ejemplo en Figura 9), información a partir de la cual se realizó el proceso de 
clasificación de cada una de las escenas. 


Figura 9 Signaturas espectrales correspondientes a los sitios de entrenamiento 
confeccionados para las imágenes del sitio 12 (Los Loros) 


Signature Comparison Chart 


suelo 

cultivo inactivo 
cultivo activo 
cultivo maduro 
veg exótica 
veg nativa 
edifi 

pavimento 
cuerpos agua 


5.9E+0001 



0 


sitio 12 band1sitio 12 band2sitio 12 bandcsitio 12 band4 final 


El proceso consiste en la clasificación de la imagen a través de un algoritmo de clasificador duro 
de tipo probabilístico denominado de “Máxima Probabilidad”, donde dependiendo de las 
características de cada escena se distribuye una función de peso que oscila entre 0 y 1 para las 
firmas utilizadas. De esta forma, se otorga una categoría a cada uno de los píxeles de la escena 
analizada, teniendo como restricción el rechazar los píxeles con menos del 95% de pertenencia 
en cada una de las clases. (Chuvieco, 2002). En términos simples, a una signatura de una 
cobertura con gran presencia dentro de una imagen (en este caso el suelo desnudo) se le 
asigna mayor “peso” probabilístico (0.5 o 50% de probabilidades de que un píxel pertenezca a 
esa cobertura de suelo), mientras que a una de escasa presencia (como los cuerpos de agua) 
se le asigna un considerablemente menor peso (en este caso cerca de 10 veces menor, con no 
más del 5% de probabilidades de ser categorizado como una cobertura). El resultado de éste 
proceso se observa en la figura 11. 
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Figura 11 Clasificación supervisada a través del criterio de máxima probabilidad para el 
sitio 12 (Los Loros) 

Máximum Likelihood Classification 
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A la imagen mostrada en la figura 3 se le realiza un filtro para eliminar píxeles que quedan 
clasificados como pequeños polígonos en polígonos de otra cobertura, generando ruidos y 
brechas temáticas. El filtro usado es el de MODA, y se utiliza un cuadrante de 3*3 píxeles. Es 
decir, si en un cuadrante de 9 píxeles, el píxel central posee un valor temático distinto al valor 
temático más repetido en ese cuadrante, se le asigna el valor más repetido. De esta manera se 
realiza una limpieza del ruido presente en la clasificación supervisada, como muestra la figura 
11 . 
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Figura 11 Imagen filtrada de la clasificación supervisada para el sitio 12 (Los Loros) 


3 x 3 Mode filter of sitio 12 max 
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Al observar la figura 9, se aprecia una considerable diferencia entre los niveles de reflectancia 
de las coberturas vegetales versus las minerales y de origen humano. Sin embargo, las 
diferencias espectrales entre cubiertas de un mismo tipo se reducen notoriamente, como en el 
caso de vegetación nativa de ribera de ríos y vegetación exótica, o entre suelos desnudos, 
construcciones y pavimento (ver figura 11), hecho que repercute en la clasificación resultante 
debió ser chequeada visualmente para minimizar la confusión entre categorías. Lo anterior se 
realizó exportando la clasificación temática filtrada a formato shape, para ser trabajada en 
ArcGIS. 


4.1 Alcances del Análisis de Teledetección en el Sector Norte 

El siguiente apartado, da cuenta de los inconvenientes o ventajas que se presentaron al 
momento de establecer las distintas clasificaciones para los sitios de estudio. 

De esta forma, el extracto presenta los alcances detectados en los sitios del 1 al 6 como 
ejemplo; no obstante, cabe mencionar que los criterios, fueron tomados para analizar toda el 
área en estudio. 

Junto con los procedimientos empleados en La clasificación espectral de las imágenes presentó 
dos enfoques, por una parte se buscó reconocer las coberturas de vegetación existentes en 
este ambiente árido, empleando como patrón los cálculos NDVI y SAVI. En este análisis para el 
caso de la vegetación de tipo agrícola se diferenció entre cultivos Maduros, en etapa de 
Crecimiento e Inactivo, este último correspondiendo a áreas que en su forma son predios 
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agrícolas pero que en términos de NDVI o SAVI, no representaban valores significativos de 
presencia de vegetación. Adicionalmente se identificaron los matorrales y pajonales, cuya 
localización esta asociada a fondos de quebradas, y la vegetación exótica presente en los 
pueblos y localidades en cada una de las escenas. En la Figura 1, se muestra el área entorno 
de la localidad de Putre, con los sitios identificados. 


Figura 12 Sitios de Entrenamiento en el Entorno de la Localidad de Putre. 



La signatura espectral promedio de las coberturas identificadas (ver Figura 13), muestra que en 
líneas generales los cultivos se diferencian de los vegetación exótica junto con matorrales y 
pajonales. La distinción más evidente es en los cultivos maduros los cuales registran la mayor 
respuesta espectral en la banda 4 (infrarrojo cercano), en tanto los cultivos en crecimiento e 
inactivos presentan una mayor reflectancias en las bandas visibles producto de la influencia del 
suelo. 
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Figura 13 Signatura Espectral en Sitios de Entrenamiento 
Localidad de Putre. 
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El segundo enfoque para establecer la presencia de vegetación cuya cobertura presenta 
diversos grados de dispersión, como es en el caso del Tolar o el Matorral desértico, fue utilizar 
como patrón de referencia la Carta de Ocupación de Tierras (COT) para el área del proyecto. Al 
respecto, es necesario indicar que dado la resolución del sensor que registra un área de 16 m2 
en el campo multiespectral, y la amplia dispersión de este tipo de vegetación en un ambiente 
árido, en términos de los algoritmos de vegetación empleados estas áreas corresponderían a 
suelos completamente desnudos. No obstante la COT, permitió tener una aproximación al área 
que cubriría cada uno de estos tipos vegetacionales, con respecto al suelo absolutamente 
desnudo. 

Lo expuesto anteriormente, se refleja en la figura 14, donde en ella se muestra una imagen 
NDVI que corresponde al sitio 1, entorno del área de Putre, en tonos verdes y azulinos se 
evidencias claramente los cultivos, en cambio los verdes claros y amarillo muestran que los 
valores de un NDVI inferior a cero domina el sector. 
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Figura 14 Resultado NDVI y COT en el Sector del Sitio 1 



De esta manera se procedió a identificar sitios de entrenamiento que permitiesen diferencias las 
áreas de suelo absolutamente desnudo de aquellas donde se presentarían coberturas 
vegetacionales asociadas a formaciones de tipo tolar o matorral desértico. En caso de Putre, la 
COT identifica un Tolar Subandino y Andino, con una cobertura de suelo desnudo de un 50 a un 
95 %. Para efectos de los sitios de entrenamiento se categorízaron como: Tolar Andino 
Disperso (75 a 90%), Tolar Subandino Disperso (90 a 95%) y Tolar Subandino Semi Disperso 
(50 a 75%). Un ejemplo en la identificación de las áreas se muestra en la figura 14. 
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Figura 15 Sitios de Entrenamiento respecto a la COT 


Andalué 


INFORME FINAL: “Levantamientos terrestres para los tramos del itinerario y paisaje cultural del 
Qhapaq Ñan-Chile, análisis y productos obtenidos de la teledetección" 26 





CONSEJO DE MONUMENTOS NACIONALES 


ANDALUE AMBIENTAL S.A. 


El análisis espectral muestra que los sectores de suelo absolutamente desnudo se diferencias 
de aquellos donde se presenta alguno de los tipos de vegetación desértica identificada. No 
obstante la diferencia entre los tipos subandino y andino son bastantes sutiles, siendo factible la 
mezcla de categorías a medida que se alejan de los sectores reconocidos por la COT. 


Figura 16 Comparación de signaturas 
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5 CARTOGRAFIA 


5.1 CONSIDERACIONES GENERALES 

El desarrollo de la cartografía de este proyecto, se desarrolló considerando la vectorización de 
los diversos elementos que fueron determinados mediante análisis satelital. Respecto de ello, 
cabe señalar algunas salvedades que generaron que finalmente algunos elementos fueran 
digitalizados directamente en pantalla: 

Las consideraciones generales de la información generada, se presenta a continuación: 

1. La unidad mínima cartográfica, con la que se debe trabajar de forma apropiada y óptima 
para información a escala 1:5.000, corresponde a 0,1 hás, esto equivale a 1000 m2 lo 
que en términos de un pixel de resolución de 16 m2 como son las imágenes kompsat 
multiespectrales, dan como resultado que al menos se necesitan 63 pixeles para definir 
una categoría, con ello sin lugar a dudas algunas clases particularmente vegetacionales 
y cuerpos de agua ínfimos no serian representados. 

2. Lo anterior, generó que se realizara un proceso de limpieza de coberturas, utilizando 
imágenes google de apoyo y fotografías del área actual tomadas durante las tareas de 
terreno. 

3. Gran parte de los impactos identificados en terreno, no presentan traducción en las 
imágenes por poseer una superficie incluso inferior a 16m 2 , por esta razón, este 
producto comprometido se entrega considerando puntos de identificación de impactos, 
los cuales fueron tomados durante las tareas de terreno. No obstante, existen unos 
pocos impactos de extensión considerable al interior de los sitios que componen el área 
de estudio que sí se observaron en las escenas y que también fueron traducidos en 
elementos puntuales. 

4. Por otra parte, los elementos lineales como vías de comunicación y cursos de agua, 
también fueron digitalizados en pantalla, puesto que por lo general (de acuerdo a lo visto 
en terreno) el ancho de los cursos no alcanza a visualizarse en las imágenes satelitales: 
bien son absorbidos por otras categorías por encontrarse sin escurrimiento superficial al 
momento de la toma de escena o, su representación en la escena es inferior a 16m 2 . 

5. Puesto que el área en estudio presenta extensas áreas con suelo predominantemente 
desnudo, es que la utilización del NVDI no permitió establecer estado de vegetación. De 
esta forma, este producto comprometido, fue solventado por la entrega de información 
referida a estado de cultivos (identificando si éstos se encontraban en crecimiento, 
maduros o si correspondían a tierras cultivadas sin uso al momento de la toma de 
escena). Para el caso de la vegetación, se utilizó la COT y la determinación de 
porcentaje de suelo desnudo que ésta presenta para así definir la densidad de la 
cobertura vegetal en cada una de las formaciones. 


5.2 PROYECCION CARTOGRAFICA 

Cada unos de los shapefiles que forman parte de los productos, se entregan en WGS 84 HUSO 
19sur. 
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5.3 ESTRUCTURA DE CARPETAS 

Cada una de las coberturas vectoriales de información se entregan en un sistema de carpetas 
donde se aprecian los siguientes apartados: 

1. Cartografía base: se presentan en esta carpeta, subcarpetas que contienen la 
información referida a vialidad, hidrografía, relieve, aldeas, toponimia y límites. 

2. Cartografía temática: Conformada por información temática como las áreas de estudio, 
hasta la información sectorizada que componen cada una de las salidas cartográficas. 
De esta forma, se presentan: los sitios arqueológicos, vegetación, áreas pobladas, 
pendientes, exposición, entre otros. 

3. Composición cartográfica: se encuentran aquí cada una de las composiciones 
cartográficas en formado MXD trabajadas con el software Arcgis9.3.1. Estas, se 
presentan en carpetas diferencidas por temáticas: zonas de entrenamiento, exposición, 
pendiente, usos de suelo, estado de cultivos y densidad de vegetación, hidrografía y 
elementos antrópicos (donde se incluye; vialidad, impactos, construcciones y vegetación 
exótica). Cabe señalar que las composiciones cartográficas en formato papel, son 
entregadas en tamaño tabloide considerando como escala de entrega de información la 
máxima probable que permite este tamaño de papel, de esta forma, las escalas de 
entrega de cada uno de los sitios en estudio, son las siguientes: 


SITIO 

ESCALA 

Sitio N 2 1 

100.000 

Sitio N 2 2 

35.000 

Sitio N 2 3 

60.000 

Sitio N 2 4 

150.000 

Sitio N 2 5 

25.000 

Sitio N 2 6 

25.000 

Sitio N 2 7 

20.000 

Sitio N 2 8 

15.000 

Sitio N 2 9 

15.000 

Sitio N 2 10 

25.000 

Sitio N 2 11 

10.000 

Sitio N 2 12 

15.000 

Sitio N 2 13 

20.000 


4. Utilidades: se adjuntan aquí las imágenes que forman las composiciones cartográficas y 
las simbologías utilizadas. 
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5.4 DICCIONARIO DE COBERTURAS 

El siguiente diccionario de coberturas, da cuenta de cada uno de los shapefiles creados como 
parte del proceso de análisis de imágenes satelitales. 

Cada una de las coberturas se presenta de forma diferenciada para cada uno de los sitios que 


forman el área de estudio. 


La información se presenta identificando 
nombre del shapefile generado: 

el nombre de los elementos que componen y el 

Nombre de elementos componentes 

Nombre de shapefile 


Vegetación nativa 

Veg_nativa. 


Comprende las formaciones de vegetación nativa presentes al interior de cada uno de los sitios 
del área de estudio. Estas fueron determinadas considerando puntos de control y Carta de 
ocupación de tierras (desarrolladas para el proyecto “Diagnóstico sobre el estado de situación 
del Qhapaq Ñan en Chile, dimensión ambiental y territorial”. 

Para el caso de la planimetría de usos de suelo, se utilizó este shp con el campo categoría. 


Actividad agrícola 


agrícola. 


Considera las tierras utilizadas para actividad agrícola, específicamente para cultivos. 

Como parte de un adicional de productos a entregar, se considera el estado en que se 
encuentran los cultivos, los cuales corresponden a: 

• Cultivos inactivos 

• Cultivos en crecimiento 

• Cultivos maduros 


Cabe señalar que el proceso de homogenización de información realizada considerando el 
conocimiento del área de estudio e imágenes google de apoyo, se concluyó que la categoría 
“suelos húmedos” que formaba parte de la entrega anterior, corresponden a tierras cultivadas 
que se encuentran en sectores cercanos a lechos de cursos de agua. 

Para el caso de la planimetría de usos de suelo, se utilizó este shp con el campo categoría. 


Construcciones 


construcciones 


Se identifican las construcciones identificadas en el área de estudio, sin categorizarlas, de esta 
forma solo presenta la categoría construcciones. 
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Suelo desnudo 

Suelo_desnudo 

Este shp se realizó considerando las clasificaciones de formación rocosa y suelo desnudo. 

Hidrografía 

cursos_agua 


cuerpos_agua 


Comprende dos coberturas de información, las cuales corresponden a cursos y cuerpos de 
agua. Cabe señalar que el desarrollo de estas coberturas se realizó considerando como base la 
información IGM disponible sobre la cual se comenzó a digitalizar en pantalla, además del 
reconocimiento de terreno y análsis visua de las imágenes KOMPSAT-2. 

Los elementos que la componen han sido categorizados de la siguiente forma: 

• Ríos 

• Quebradas 

• Cuerpos de agua 

Cabe señalar que para el caso del sitio 10, no se reconocieron cursos de agua utilizando 
imágenes google, no existían puntos de control (puesto que el elemento no fue identificado) y 
no se reconocieron mediante análisis de imágenes satelitales, por lo que la cartografía asociada 
para este sitio, corresponde a la visualizacion de riadas que se reconocieron mediante filtros 
lineales en Idrisi. 


Vialidad 


vialidad 


Comprende la identificación de la vialidad que se presenta en los sitios de entrenamiento, éstas 
fueron desarrolladas considerando como base la información OGM disponible y reconocimiento 
de terreno. 


Como adicional a este producto, se presenta la categorización de acuerdo a tipo de carpeta, los 
cuales corresponden a: 

• Pavimento 


• Tierra 


Impactos 


impactos 


La información contenida en este shapefile corresponde a los impactos identificados en terreno 
durante la campaña de terreno y que corresponden a los sitios de entrenamiento. La evidencia 
de impactos de mayores dimensiones, fue determinada mediante el análisis de las imágenes 
satelitales. 
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Areas pobladas 

Area_poblada 


Corresponde a polígonos que delimitan las áreas pobladas de acuerdo a la densidad de las 
construcciones. 

Dado el tamaño de información para el sitio 4, no fue posible cruzar la clasificación con las 
áreas pobladas, por esta razón, se adjunta el shp de áreas pobladas. 


Areas industriales 

Areajndustrial 




Corresponde a aquellas áreas en donde la actividad industrial se desarrolla actualmente. 


Pendiente 


pendienteclas 


El desarrollo de esta información se llevó a cano considerando las distintas escenas analizadas. 


La clasificación que se utilizó para cartografiar y diferenciar corresponde a los criterios 
establecidos en el Manual Técnico N e 23 de CONAF “Método para la Planificación del Manejo de 
Areas Protegidas”, estas categorías corresponden a: 


Calificación 

MUY ADECUADA 
ADECUADA 

MEDIANAMENTE ADECUADA 
POCO ADECUADA 
NO ADECUADA 


Descriptor 

Pendientes entre 0 y 15 2 
Pendientes entre 15 y 30 2 
Pendientes entre 30 y 45 2 
Pendientes entre 45 y 60 2 
Pendientes superiores a 60 2 


Exposición 


exposicionclas 


Este shapefile de información, se desarrollo considerando la exposición en forma cardinal, 
siendo las categorías de ésta, las siguientes: 

• Norte 

• Noreste 

• Este 

• Sureste 

• Sur 

• Suroeste 

• Oeste 

• Noroeste 

• Plano 
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Zonas de entrenamiento 


zentrenamiento 


Considera los distintos puntos utilizados en el proceso de análisis satelital, estos abarcan a 
aquellos puntos tomados en terreno, como aquellos que fueron identificados en google earth. 


Con la información proveniente de “agrícola” y “vegetación nativa” se compusieron los planos de 
estado y densidad de vegetación, utilizando para ello el campo clasif para distinguir el estado de 
los cultivos y la densidad de la vegetación. 
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5.5 DICCIONARIO DE CAMPOS 

Cada uno de los shapefile generados, presenta los siguientes campos: 

• Tipo: corresponde a la identificación de cada uno de los registros, el cual puede 
corresponder por ejemplo a pendiente no adecuada, Pajonal muy denso, 
construcciones, entre otros. 

• Clasif: campo donde se identifica una primera agrupación de información de acuerdo a 
su utilidad para cartografiarse, por ejemplo, los registros identificados como: pajonal 
poco denso y matorral ripario. 

• Categoría: identifica la asociación de diversas clasificaciones, por ejemplo: la categoría 
suelo desnudo, asocia a varias clasificaciones como son lecho de río, caminos 
pavimentados detectados por procesamiento de imágenes satelitales, vegtaciones cuyo 
porcentaje de suelo desnudo supera el 90% y pavimento desértico. 

• Area_ha: identifica el área en hectáreas para cada uno de los registros tipo polígonos 

• Longm: identifica la longitud de cada uno de los elementos tipo línea que conforman los 
distintos shp. 

• Norte/este: utilizado para los shp de elementos puntuales, donde se identifica la 
coordenada norte y este que presentan. 
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5.6 LISTADO DE PLANOS 


PLANO N° 

TEMATICA 

AREAS DE ENTRENAMIENTO 

1 

Areas entrenamiento sitio N 2 1 

2 

Areas entrenamiento sitio N 2 2 

3 

Areas entrenamiento sitio N 2 3 

4 

Areas entrenamiento sitio N 2 4 

5 

Areas entrenamiento sitio N 2 5 

6 

Areas entrenamiento sitio N 2 6 

7 

Areas entrenamiento sitio N 2 7 

8 

Areas entrenamiento sitio N 2 8 

9 

Areas entrenamiento sitio N 2 9 

10 

Areas entrenamiento sitio N 2 10 

11 

Areas entrenamiento sitio N 2 11 

12 

Areas entrenamiento sitio N 2 12 

13 

Areas entrenamiento sitio N 2 13 

PENDIENTES 

14 

Clasificación Pendientes sitio N 2 1 

15 

Clasificación Pendientes sitio N 2 2 

16 

Clasificación Pendientes sitio N 2 3 

17 

Clasificación Pendientes sitio N 2 4 

18 

Clasificación Pendientes sitio N 2 5 

19 

Clasificación Pendientes sitio N 2 6 

20 

Clasificación Pendientes sitio N 2 7 

21 

Clasificación Pendientes sitio N 2 8 

22 

Clasificación Pendientes sitio N 2 9 

23 

Clasificación Pendientes sitio N 2 10 

24 

Clasificación Pendientes sitio N 2 11 

25 

Clasificación Pendientes sitio N 2 12 

26 

Clasificación Pendientes sitio N 2 13 

EXPOSICION 

27 

Exposición Sitio N 2 1 

28 

Exposición Sitio N 2 2 

29 

Exposición Sitio N 2 3 

30 

Exposición Sitio N 2 4 

31 

Exposición Sitio N 2 5 

32 

Exposición Sitio N 2 6 

33 

Exposición Sitio N 2 7 

34 

Exposición Sitio N 2 8 

35 

Exposición Sitio N 2 9 

36 

Exposición Sitio N 2 10 

37 

Exposición Sitio N 2 11 
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38 

Exposición Sitio N 2 12 

39 

Exposición Sitio N 2 13 

USOS DE SUELO 

40 

Uso de suelo Sitio N 2 1 

41 

Uso de suelo Sitio N 2 2 

42 

Uso de suelo Sitio N 2 3 

43 

Uso de suelo Sitio N 2 4 

44 

Uso de suelo Sitio N 2 5 

45 

Uso de suelo Sitio N 2 6 

46 

Uso de suelo Sitio N 2 7 

47 

Uso de suelo Sitio N 2 8 

48 

Uso de suelo Sitio N 2 9 

49 

Uso de suelo Sitio N 2 10 

50 

Uso de suelo Sitio N 2 11 

51 

Uso de suelo Sitio N 2 12 

52 

Uso de suelo Sitio N 2 13 

ESTADO CULTIVOS Y DENSIDAD VEGETACION 

53 

Estado de cultivos y densidad de vegetación Sitio N 2 1 

54 

Estado de cultivos y densidad de vegetación Sitio N 2 2 

55 

Estado de cultivos y densidad de vegetación Sitio N 2 3 

56 

Estado de cultivos y densidad de vegetación Sitio N 2 4 

57 

Estado de cultivos y densidad de vegetación Sitio N 2 5 

58 

Estado de cultivos y densidad de vegetación Sitio N 2 6 

59 

Estado de cultivos y densidad de vegetación Sitio N 2 7 

60 

Estado de cultivos y densidad de vegetación Sitio N 2 8 

61 

Estado de cultivos y densidad de vegetación Sitio N 2 9 

62 

Estado de cultivos y densidad de vegetación Sitio N 2 10 

63 

Estado de cultivos y densidad de vegetación Sitio N 2 11 

64 

Estado de cultivos y densidad de vegetación Sitio N 2 12 

65 

Estado de cultivos y densidad de vegetación Sitio N 2 13 

HIDROGRAFIA 

66 

Hidrografía Sitio N 2 1 

67 

Hidrografía Sitio N 2 2 

68 

Hidrografía Sitio N 2 3 

69 

Hidrografía Sitio N 2 4 

70 

Hidrografía Sitio N 2 5 

71 

Hidrografía Sitio N 2 6 

72 

Hidrografía Sitio N 2 7 

73 

Hidrografía Sitio N 2 8 

74 

Hidrografía Sitio N 2 9 

75 

Hidrografía Sitio N 2 10 

76 

Hidrografía Sitio N 2 11 

77 

Hidrografía Sitio N 2 12 

78 

Hidrografía Sitio N 2 13 

ELEMENTOS ANTROPICOS 

79 

Elementos antrópicos sitio N 2 1 
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80 

Elementos antrópicos sitio N 2 2 

81 

Elementos antrópicos sitio N 2 3 

82 

Elementos antrópicos sitio N 2 4 

83 

Elementos antrópicos sitio N 2 5 

84 

Elementos antrópicos sitio N 2 6 

85 

Elementos antrópicos sitio N 2 7 

86 

Elementos antrópicos sitio N 2 8 

87 

Elementos antrópicos sitio N 2 9 

88 

Elementos antrópicos sitio N 2 10 

89 

Elementos antrópicos sitio N 2 11 

90 

Elementos antrópicos sitio N 2 12 

91 

Elementos antrópicos sitio N 2 13 
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